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Abstract 

Several topics from the wide field of QCD studies in Deep-inelastic 
ep Scattering at HERA are addressed. They include QCD analyses of 
the inclusive cross section with the determination of ag and the proton 
gluon density from the F2 scaling violations, and the determination of 
the longitudinal structure function Fl. QCD analyses of inclusive jet 
and dijet data are also presented. Finally jet substructure and three-jet 
production are discussed. 


1 Introduction 

Studies of Deep-inelastic lepton-nucleon scattering (DIS) have played a fun¬ 
damental role in establishing QCD as the strong interaction theory and in 
exploring the parton structure of the nucleon. With the advent of HERA, in 
which electrons or positrons of 27.5 GeV energy collide with protons of 820 
GeV (in the last years 920 GeV), the tests of QGD have been extended by 
several orders of magnitude with respect to the range in Bjorken-x and in Q^, 
the squared momentum transfer between lepton and nucleon. The early hxed 
target experiments observed scaling violations, i.e. the variation of the struc¬ 
ture functions with Q^. The scaling violations are well described by QGD, in 
which they are related to the gluon density in the proton, and to the strong 
interaction coupling constant, 

At HERA the low x region was experimentally explored for the first time, and 
the first measurements [jl[] of the proton structure function F 2 {x,Q‘^) revealed 
a steep rise of F 2 towards low x values, indicating a high gluon density in the 
proton at low x. A key question is then the validity of the DGLAP evolution 
equations]^] at low x values, since in the DGLAP evolution higher order terms 
proportional to ag ■ ln{l/x) are neglected. One expects that at some value 
of X non-linear gluon interaction effects will become important, damping the 

*Talk given at the Intern. Conference on New Trends in High Energy Physics, 
22-29.9.2001, Yalta, Crimea, Ukraine. 
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rise of the cross section in accordance with unitarity requirements. While this 
question cannot yet be answered, the low x region remains a key area for future 
QCD studies at HERA. 



Figure 1: Diagram for 
lowest order ep DIS 


At the same time as precision measurements of the 
inclusive DIS cross section and the proton struc¬ 
ture functions are being pursued, QCD tests of 
processes of higher order in are also being in¬ 
vestigated. These tests involve the study of jet 
production in DIS, a held for which HERA is well 
suited with its large ep CMS energy a/s = 300—320 
GeV. This experimental work is intimately con¬ 
nected with theory development, since tests of 
higher order jet production need the corresponding 
processes to be quantitatively estimated. 


2 The Inclusive DIS Cross Section and the 
Proton Structure 

Inclusive ep neutral current (NC) DIS is to lowest order, described by 

the diagram in Fig. 1, in which a virtual boson is emitted by the electron and 
scatters off a parton in the proton. In the one-photon exchange approximation, 
the double differential cross section can be written as 

^ - y^Fx(x,Q^)), = 1 + (1 - vf. 

The inelasticity y = jxs represents the energy fraction transferred to the 
proton in the scattering process. and are the (unpolarized, electro¬ 
magnetic) proton structure functions and contain the information about the 
momentum distribution of partons in the proton. They are not calculable in 
theory and have to be measured. However, when measured at a given value, 
their evolution in can be predicted in theory. 

The NC DIS cross section can also be written as 

= V{y){aT F e[y)aL) 

where the ep scattering process is considered as the interaction of a flux of 
virtual photons with the proton. Here r(j/) = Y^aji^'xQ^x') is the flux factor 
and e{y) = 2(1—|/)/K|_ dehnes the photon polarization, ax and ai are the cross 
sections of the interaction of transverse and longitudinally polarized photons 
with the proton. These cross sections are related to the structure functions: 

F 2 {x, Q^) = ^{arix, Q^) + aiix, Q^)) and Fl(x, QQ = QQ. 

Due to cross section positivity, the relation 0 < Fl < F 2 is obeyed. 

In the QPM world, without gluons, ai = 0, since longitudinally polarized 
photons do not interact with massless spin 1/2 partons. Thus, in QPM also 
Fl = 0[^. In QCD quarks radiate gluons and interact through gluon exchange. 
Radiated gluons in turn can split into quark-antiquark pairs (“sea quarks”) or 
gluons. The gluon radiation results in a transverse momentum component 
of the quarks, which can now also couple to longitudinally polarized photons. 
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Thus, ai and Fl get non-zero values. Due to its origin, Fl is directly dependent 
on the gluon distribution in the proton and therefore the measurement of Fl 
provides a sensitive test of QCD at low x values. In fact, the low x region 
cannot be understood without measuring Fl precisely. 

The ep cross section is usually written 
in the “reduced” form, in which the 
dependence due to the photon propa¬ 
gator is removed, 

ar = F2{x,Q^)-^FL{x,Q^). 
Since the contribution of the longitu¬ 
dinal structure function Fl to the cross 
section can be sizeable only at large val¬ 
ues of y, in a large kinematic range the 
relation cr^ ~ F 2 holds to a very good 
approximation. 

The HI and ZEUS collaborations 
have recently presented high statistics 
measurements 0, ||, || of the NC DIS 
cross section and extracted F 2 in the 
kinematic range 1.5 — 2.7 < < 

30000 GeV^ and 3-6-10-5 < a; < 0.65. 
The ZEUS data are shown in Fig. 2. 
The strong (positive) scaling violations 
at low X values, due to the increase of 
the gluon density ~ cusXg at 

low x), are clearly seen. At large x val¬ 
ues negative scaling violations appear 
(alfV ^ “"^2 at large x). 

The ZEUS and HI data are in very good agreement, both with each other 
and, at the largest x values where data overlap, with the earlier hxed target 
experiments. 

Both collaborations have subjected their F 2 data to extensive QCD analyses in 
next-to-leading-order (NLO), extracting the parton density functions (PDFs) 
and OsP, 0]. While the general analysis and £t strategy is similar in both 
collaborations, there are also many differences in the details, e.g. in the density 
parametrizations and in the treatment of flavour number and in the use of the 
hxed target data. An extensive discussion can be found in P). 

Fig. 2 shows that the NLO QCD hts, which are based on the DGLAP evolution, 
describe the data very well over no less than four orders of magnitude and 
down to surprisingly low of a few GeV^. The gluon density extracted in 
the NLO QCD hts is shown for both experiments in Fig. 3 for three values 
of Q^. Within the error bands there is reasonable agreement. Like T 2 , the 
gluon density increases towards low x values, and the increase is steeper with 
increasing Q^. 


ZEUS 



Figure 2: F 2 vs Q"^, for hxed val¬ 
ues of X. Data from ZEUS and 
hxed target experiments. The 
curves show the NLO QCD ht to 
the data. 
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Figure 3: Gluon density for three 
different values of Q^. The error 
bands show the experimental and 
(for HI) also the total (including 
theoretical) uncertainty. 
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Figure 4: as values obtained in 
the NLO QCD analyses of F 2 
data. Also shown are as val¬ 
ues obtained in a recent NNLO 
analysis]^ as well as other 
values obtained at NNLO level. 


The as values obtained in the NLO hts of HI and ZEUS are shown in Fig. 4[|. 
In the HI case also the error due only to the experimental uncertainty is shown. 
Thus, the theoretical uncertainty dominates the total error. The major part of 
this error is due to the choice of renormalization scale, and to a lesser extent, 
also to the choice of the factorization scale. This is expected to change once 
next-to-next-to-leading-order (NNLO) QCD calculations are available. Fig. 4 
also shows the result of a recent NNLO analysis of ep DIS data[]^, in which 
the moments of F 2 are htted to obtain an impressingly small error on 
Comparing the ep DIS as result with the result obtained when using the more 
imprecise uN scattering data, or with the as values obtained in NNLO for 
r-decays, T-decays, or in analysis of F(Z), it is obvious that the ep DIS data 
are very competitive in precision determinations of as- Note that the analysis 
in used earlier data on ^ 2 ( 0 , and not yet the most recent, very precise 
data. 

Both collaborations are working on further reduction of the experimental errors 
on the F 2 data, errors which are completely dominated by systematics, in 
particular the uncertainty on the calorimetric energy scale. Clearly, theoretical 
progress is essential in order to proht from future, still better experimental 
precision in the F 2 data. 


3 The Longitudinal Structure Function 

A fully independent measurement of Fl at HERA, in a wide range of x and 
is only possible via a substantial variation of the beam energies; at any 

^Figure courtesy M. Erdmann, LP2001, Rome 2001 Q. 
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given values of x and Q‘^ the difference of the reduced cross section, measured 
at two values of the CMS energy s = Q'^/xy, is a direct measure of 

- (Tr,2 = [(2/V^+)l - (l/V^+)2] ' Fl- 



Figure 5: Reduced cross 
section vs. x, for = 
12 GeV^. The curves 
are due to the HI NLO 
QCD £t and represent 
F 2 and the cross section. 


This measurement is not yet performed but is 
part of the HERA H physics program]^. While 
waiting for this to happen, the HI collaboration 
has performed several indirect extractions of F^. 
The methods of extraction are based on the di¬ 
rect measurement of the cross section in the high 
y region, and on the extrapolation of the precise 
knowledge of F 2 at low values of y, into the re¬ 
gion of high y. Cross section measurements of 
the NC DIS process in the high y region pose a 
challenge to the HERA experiments, since low 
energy (down to 3 GeV is achieved) scattered 
electrons have to be triggered, identihed and well 
measured. The detailed understanding of the 
detector and of the photoproduction background 
is essential and an important part of the mea¬ 
surement is therefore the improvement of both 
calorimetry and the detection of charged tracks 
in the HI central and backward regions. The 
measurement of track charge and momentum at 
all scattering angles is particularly important in 


estimating the amount and shape of the background from photon conversions 
(with photons from the 7 r° —> 77 decay), through the identihcation of ’’wrong” 
charge electrons. 

Two methods are used for the extraction of Fl from the cross section mea¬ 
surements, namely the ’’subtraction” method applied at > 10 GeV^, and 
the ’’derivative” method, applied for < 10 GeV^. The subtraction method 
is illustrated in Fig. 5, showing the reduced cross section as function of x at 
hxed = 12 GeV^. At the lowest x values (corresponding to high y values) 
the cross section falls below the F 2 curve, extrapolated via the HI NLO QGD 
£t. Thus, Fl is obtained from the difference F 2 — ar- This method was hrst 
explored in |T^ and in the recent analyses, using data from 1996-97^ and 
1998-2000 |T^|, Fl is extracted up to values of 700 GeV^. Fig. 6 shows the 
dependence of Fl, separately for e~^p and e~p data. As it should be, Fl is 
independent of the lepton beam. One also sees that Fl is clearly different from 
the extreme possibilities, F^ = 0 or F^ = F 2 . Within the precent precision 
of the data there is good agreement with the QGD expectation, as projected 
using the HI NLO fit. 

In the second method of extracting Fl, the derivative of the reduced cross 
section with respect to In y is formed: 



dPL 
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Figure 6: vs. (or, equiv¬ 

alently, vs. X, upper scale), for 
fixed y = 0.75, for e^p and e~p 
data. The curves show the HI 
NLO QCD £t, obtained for sev¬ 
eral assumptions on F^. The 
shaded band shows the expecta¬ 
tion for Fl and its uncertainty 
from the QCD £t. 


For y ^ 1 the cross section derivative 
tends to the limit ( 9 ^ 2/9 In j/)q 2 — 
2 ■ Fl, neglecting the contribution 
from the derivative of Fl. At largest 
y the Fl contribution dominates the 
derivative of the reduced cross section 
ar- This is in contrast to the influence 
of F7 on ar which is dominated by the 
contribution of F 2 for all y. A further 
advantage of the derivative method is 
that it can be applied down to very 
low ~ 1 GeV^ where a QCD de¬ 
scription of F 2 {x,Q^) is complicated 
due to higher order and possible non- 
perturbative corrections. The mea¬ 
sured cross section derivative is shown 
in Fig. 7 as a function of y, in several 
bins of Q^. Since for a given Q'^ value, 
F 2 ~ x~^ ~ and since A is small 

for low values of is linear in 

\ny to good approximation. This is 
clearly seen in Fig. 7. In each bin 
straight line hts were made to the de¬ 


rivative data for y < 0.3. The line hts describe the data very well and the 
extrapolation of the straight line was taken to represent the contribution of F 2 
at high y. The small contribution of dFL/d\\iy to the derivative was corrected 
for by using NLO QCD and the correction was included in the overall error of 
the measured Fl- The derivative method was used in the analysis of the low 
data, both in the 1996-97 data|^ and in the 1999-2000 datall^. Since the 
latter data have higher p beam energy, 920 Gev as compared to 820 GeV in the 
earlier data, the accessible x range could be extended to even smaller values. 
It is also worthwhile to note that the HI NLO fit is based on the 1996-97 data; 
when using this £t in the analyses of the later data proper account was taken 
in order to correct for the difference in p beam energy. 
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Figure 7: Measurement of 
{dar/d\n.y)Q 2 vs. y, in 
bins of Q^. The curves 
represent the NLO QGD 
£t to the HI 1996/97 data 
for y < 0.35 and > 3.5 
GeV^, calculated for sev¬ 
eral assumptions on Fl- 


6 




















a 



-4 -2 

10 10 1 


c 

g 

o 

SI 

"o 

O 


X 


X 


Figure 8: Longitudinal Structure Function Fl vs. x and in bins of 
Q^, obtained by HI and fixed target experiments. Error bands are 
due to experimental (inner) and model (outer) uncertainties using 
the HI NLO QCD fit to the HI 1996/97 data for y < 0.35 and 
> 3.5 GeV2. 


Fig. 8 gives an overview of the current HI data on Fl{x, Q'^) in the range 
2.2 — 700 GeV^I^ 1^. The data extend the knowledge of Fl into the region of 
low X, much beyond the range of the hxed target experiments. The increase of 
Fl{x, g^) towards low x is consistent with the NLO QGD calculation, reflecting 
the rise of the gluon density in this region. The values of Fl{x,Q‘^) are thus 
severely constrained by the present data, unless there are deviations from the 
assumed extrapolation of F 2 into the region of large y corresponding to the 
smallest x. 

Fl has also been calculated in the ZEUS NLO QGD fit[0; the calculation is 
consistent with the HI data. 
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4 Inclusive Jets and Dijets in DIS 


The fully inclusive DIS process, shown in Fig. 1, is 
of zero order in There is direct sensitivity only 
to the quark PDF’s, and the determination of as and 
gluon density is indirect, via the scaling violations of 
F 2 , measured over a large range of x and Q^. 

The DIS processes shown in Figs. 9a and b, QCD 
Compton (QCDC) scattering and Boson Gluon Fu¬ 
sion (BGF), are of order Ug and due to the latter 
contribution there is direct sensitivity to the gluon 
density of the proton. At high enough energy, the in¬ 
volved final state partons manifest as jets of hadrons. 
The multi-jet final state can be characterized by the 
jet mass, Mjj, and the variable ^ = XBj{^ + 

The dijet mass Mjj gives the CM energy of the boson- 
parton reaction and the fractional momentum x (the 
longitudinal momentum fraction of the proton car¬ 
ried by the parton specified by the PDFs) is given by 
1 ^, which is much larger than xsj if Mjj is large. 

Two hard scales enter in these diagrams, and 
Erjet- Studies of the dynamics of multi-jet pro¬ 
duction are usually performed in the Breit frame, 
where the virtual boson interacts head-on with the 
proton]^. In lowest order, the parton from 

the proton is backscattered into opposite direction, and no transverse energy 
is produced. Thus, appearance of jets with large Et can only be explained by 
hard QCD processes of at least 0{as), and E^jet is then the physical scale at 
which e.g. hard gluon radiation is resolved. Experimentally, jets are usually 
identified using a fcr algor it hm|p.ti||. 

The HI and ZEUS collaborations have recently presented high statistics studies 
of inclusive jet and dijet production in NC DIS[^, |^, The data span 

the kinematic range 5 — 10 < < 10 — 15 ■ 10^ GeV^ and 5 — 7 < ETjet < 

60 GeV. NLO QGD calculations in general describe the data very well, over 
almost the whole range of and E-rjet- This is exemplihed in Fig. 10, 
showing the HI inclusive jet cross section as function of E-rjet and Q^, and in 
Fig. 11, showing the ZEUS dijet cross section as function of Q^. The largest 
uncertainty in the QGD calculations stems from the uncertainty in the choice 
of renormalization scale i^r, which is taken either as Q or as mean jet transverse 
energy, Erjet- The uncertainty is largest for low values (Fig. 11b), and cuts 
are made at 150 (HI) and 470 (ZEUS) GeV^ in the QGD analyses of the data. 
At large values of the experimental uncertainties are also smaller, and this 
is true as well for the hadronic (parton-to-hadron) corrections to the NLO 
calculations, as shown in Figs. 10 and 11c. 

Once the good agreement of the QGD calculation with the jet data has been 
established, one can proceed to perform a QGD analysis, determining and 




Figure 9: Diagrams 
for QGD Gompton 
and Boson Gluon 
Fusion processes. 
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the gluon density. The analysis strategy is as follows: 

1. Fit the jet data to the NLO QCD predictions, using PDFs obtained 
from global, external £ts[^, With the PDFs externally hxed, there 
is only one free parameter, namely as- Note that the PDFs themselves 
depend implicitly on and that this complication was properly taken 
into account in the hts. 

2. After establishing that the htted as value agrees well with other, external 
measurements, is hxed (e.g. to the world average) and the jet data are 
then htted in order to extract the PDFs, in particular the gluon density 
of the proton. 

3. Finally, a global, simultaneous ht of both as and the PDFs can be per¬ 
formed. This would be a more independent test of QCD with the data. 
However, although in this global ht the quark PDFs emerge as identical 
to those resulting from the ht in step (2), the simultaneous ht fails to 
produce meaningful results for as and the gluon density. The reason is 



Figure 10: Inclusive jet cross sec¬ 
tion vs. in several inter¬ 

vals of Q^. Also shown are NLO 
QCD calculations, with and with¬ 
out hadronic corrections. 


ZEUS 



Figure 11: Dijet cross section 
vs. Q^. The shaded band repre¬ 
sents the calorimeter energy scale 
uncertainty. The curves show 
NLO QCD calculations for sev¬ 
eral choices of and proton PDF 
sets, b) Ratio data/theory show¬ 
ing in addition the ehect of the 
theoretical scale uncertainty, 
c) Parton-to-hadron correction. 


the strong anti-correlation between as and the gluon density, which can 
be understood from the fact that in the phase space region considered, 
jet production is dominated by the gluon contribution, which enters in 
the cross section as the product ag ■ xg{x). 
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Figure 12: a* running in 2 scales, Erjet and Q. Uncertainties due to 
calorimeter energy scale and to theory are indicated in the ZEUS plot. 
Upper curves show the RGE prediction. 
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Figure 13: Gluon density in the 
proton, determined in a combined 
QGD £t using inclusive DIS , in¬ 
clusive jet and dijet cross sections. 
The error band includes all ex¬ 
perimental and theoretical uncer¬ 
tainties, also that of as{Mz)- 


Both HI and ZEUS have performed 
step 1 in this strategy, and the resulting 
as measurements are shown in Fig. 11. 
Note that the “running” of is here 
clearly seen in both scales, Q and ETjet, 
and the running is moreover seen within 
each single experiment. The running 
is consistent with the renormalization 
group equation (RGE). The comparison 
of the measurements with other HERA 
measurements, and with the world av¬ 
erage values is given in Fig. 15. The 
fact that the ag values obtained from 
the DIS jet data agree very well with 
other measurements, in particular with 
measurements from processes which do 
not involve hadrons in the initial state, 
like annihilation to hadrons, can be 
taken as proof of the validity of pertur¬ 
bative QGD at NLO in DIS jet produc¬ 
tion. 


Step 2 in the strategy was performed by HI, and the gluon density extracted 
from the data is shown in Fig. 13. The fit includes the combined cross sec¬ 
tions of inclusive jet data, dijet data and inclusive DIS. The result is in good 
agreement with recent global analyses pT|, |2^ , p3| . For a detailed comparison 
of this result with the gluon densities obtained from the E 2 scaling violations, 
see Pl . 
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5 Jet Substructure 



Figure 14: a) Integrated jet shape and b) mean subjet multiplicity as 
functions of ETjet- NLO QCD calculations are shown for three different 
values of ag- 


At sufficiently high transverse jet energy, -Ft jet (or fragmentation ef¬ 

fects become negligible and both the shape and other features of the internal 
jet structure are expected to be calculable in perturbative QCD. The ZEUS 
collaboration has recently presented |]25| several results based on studies of jet 
substructure, using the variables p6|, ^ “integrated jet shape” \l/(r) and “mean 
subjet multiplicity” {usubjet)- 


'h(r) 


1 

Njet Ei^\r=R) 


and sub jet) 


ll „ . 

N. Z^l itsubjet- 


The radius r is dehned in (j) — t] space, where the jet search is performed. The 
subjets within a given jet are found by repeating the jet algorithm with smaller 
resolution scale. 

Among the results presented in ||^ are 


• The average subjet multiplicity in jets in Charged Current (CC) and 
NC events is found to be similar. Since the jets in CC and NC DIS are 
predominatly quark initiated, the similarity in jet substructure indicates 
that the pattern of parton radiation within a quark jet is independent of 
the specihc hard scattering process. 


• In a dijet event sample, c-quark jets were tagged through the identihca- 
tion of a D* meson. The internal structure of the charm induced jets is 
found to be similar to that of the quark jets in NC DIS. Since the latter 
are dominantly light quark initiated, one can conclude that the evolu¬ 
tion of the outgoing partons, which determines the internal structure of 
the jet, is independent of the hard subprocess from which the outgoing 
partons originate. 
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• Using the internal jet structure of the tagged, charm induced jets in the 
dijet event sample, and comparing with the internal jet structure of the 
total dijet sample, it was possible to extract the internal jet structure 
of gluon jets. The prediction of QCD, that gluon jets are broader, and 
contain more subjets, is nicely conhrmed. 


HERA Measurements 
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(J. Phys. G26 (2000) R27) 
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Figure 15: Summary of ag mea¬ 
surements at HERA (given at the 
Z mass), using inclusive DIS data 
and NLO QCD hts, inclusive jets 
and dijets, and jet shape and sub¬ 
jet multiplicity. Also shown are 
the world average values (PDG 
2000, S. Bethke). 


The internal jet structure is sensitive 
to ttg beyond leading order. This sensi¬ 
tivity is demonstrated in Fig. 14 a and 
b, where the integrated jet shape (for 
r = 0.5) and mean sub jet multiplicity 
is shown as functions of Etjci- Note 
that the jets become narrower, and that 
the mean subjet multiplicity decreases 
as ETjet increases. The NLO QCD 
calculation describes both jetshape and 
mean subjet multiplicity well, but the 
calculation also varies strongly with dif¬ 
ferent values of ag. It is thus clear that 
the jet substructure data can be used 
to determine ag. The method is simi¬ 
lar to the one used in the QCD analysis 
of the inclusive and dijet data described 
above. External sets of PDFs are used, 
and the implicit dependence on ag in 
these sets is properly taken account of 
in the htting. 

The results of the ag determination 
using the jet substructure are shown 
in Fig. 15, together with other ag 
measurements from HERA. These new 
measurements agree well with the oth¬ 
ers, which is proof of the consistency of 
the NLO QCD calculations also where 
the internal jet structure is concerned. 


A general remark about the measurements in Fig. 15 can be made: The exper¬ 
imental errors in these measurements are small and comparable to the error on 
the world averages PB|, The total errors on the ag measurements at HERA 
are in fact everywhere dominated by the theoretical uncertainty. As already 
noted above, this situation is expected to change when NNLO calculations 
become available. 
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The Three-Jet Final State in DIS 




Figure 16: a) Inclusive three-jet 
cross section vs. Q^. Also shown 
is the ratio of data to theoretical 
prediction, b) Ratio of inclusive 
three-jet and dijet cross sections 
vs. Q^. LO and NLO QCD cal¬ 
culations are also shown, includ¬ 
ing the effects of variation in the 
latter of as{Mz), renormalization 
scale fir and proton gluon density. 


While the inclusive jet and dijet cross 
sections are directly sensitive to QCD 
effects of order 0{as), the three-jet 
cross section in DIS is already pro¬ 
portional to oPg in leading order in 
perturbative QCD. This higher sen¬ 
sitivity to as and the greater number 
of degrees of freedom of the three-jet 
hnal state allow the QCD predictions 
to be tested in more detail in three-jet 
production. The HI collaboration has 
recently presented a study of DIS 
three-jet events, covering the kine¬ 
matic range 5 < < 5000 GeV^ and 

three-jet masses 25 < M^j^t < 140 
GeV. The dependence of the cross 
section is shown in Fig. 16a, together 
with QGD LO and NLO calculations, 
the latter with and without hadronic 
corrections. The NLO calculation, 
which is due to the recently available 
program NLOJET[^], describes the 
data well over the whole kinematic 
range, as is also seen in the ratio of 
data to theory. In the latter plot are 
also shown the theory uncertainty due 
to variation of the gluon density, the 
renormalization scale fir = and 
as- As seen, at large > 50 GeV^ 
the as variation gives the largest un¬ 
certainty. 

Since both dijet and three-jet pro¬ 
duction is dominated by gluon in¬ 
duced processes, the uncertainty of 
the gluon density can be minimized 
by taking the ratio R 3/2 of three-jet 


and dijet cross sections, at the same values of x and Mn-jet- R can be shown, 
using the QGD calculations, that dijet and three-jet production involves the 
same gluon fraction, at similar values. Furthermore, many experimental 
and systematic errors cancel in the ratio. As is evident in Fig. 16b, R 3 / 2 , which 
is directly proportional to as, is experimentally measured and theoretically 
calculated with small uncertainties over the whole range. Thus, given 
better statistics, a very sensitive test of QGD will be possible with the three- 
jet data, including a precision measurement of as- 
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Figure 17: Normalized distributions of cos^s 
and angle "03 in the three-jet CMS at low 
and high Q^. Solid and dashed curves show 
NLO and LO QCD calculations, dotted curves 
show a three-jet phase space model. 


The topology of the three- 
jet hnal state offers a test 
of QCD through the an¬ 
gular distributions of 6*3 
and V’ 3 ) where 63 is the 
angle between the highest 
energy jet and the proton 
beam, and -03 is the angle 
between the two planes 
formed by the highest en¬ 
ergy jet and the pro¬ 
ton beam, and by the 
three jets, respectively. 
These angles are shown 
in Fig. 17, for two 
ranges. Both LO and 
NLO QCD calculations 
describe these normalized 
distributions well, while 
a phase space model fails 
the description. The dis¬ 
tributions show that the 
jets tend to be aligned 
with either the photon or 


the proton, i.e. the Bremsstrahlung nature of the process (coherence property 
of QCD) is conhrmed. Similar distributions were previously also observed by 
the ZEUS collaboration in three-jet photoproduction [I 


Conclusions 

The DIS data are well described by NLO QCD and the DGLAP evolution, 
over a large range of and Bjorken-x. This is true for the inclusive DIS cross 
section, with the structure functions F 2 and F^, as well as for the inclusive 
jet production and the exclusive dijet and three-jet production and for the 
internal structure of jets. 

In many areas the data have reached a high experimental precision, and 
progress in the tests of perturbative QCD depends crucially on further progress 
in the theory, where the NNLO calculations for the DIS processes are awaited. 
This is particularly true for the precision determination of ag at HERA, using 
DIS data, where the experimental precision is already at level with the error 
on the world average, and is expected to improve even further. 

The running of ag in accordance with the RGE is seen within each single 
experiment, and in both hard scales, Q and E^jet- 

Some areas, like jets at highest and Ft, jet, and the three-jet hnal state, 
are still statistically limited. The HERA H running period will bring a huge 
improvement of statistics in the coming years. However, in the analyses de- 


14 
















scribed in this report, the experiments have in many cases so far only used the 
data taken in 1995-97, which constitute only about 1/3 of the total HERA I 
data. The remaining 2/3 of the HERA I data, from the years 1998-2000, is 
currently being recalibrated and reprocessed, and will soon be available. 
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